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Abstract 

A new class of soliton-like solutions is derived for the Grad-Shafranov (GS) equations. A math- 
ematical analogy between the GS equation for MHD equilibria and the cubic Schrodinger (CS) 
equation for non-linear wave propagation forms the basis to derive the new class of solutions. The 



. soliton-like solutions are considered for their possible relevance to astrophysics and solar physics 
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t^J- problems. We discuss how a soliton-like solution can be generated by a repetitive process of mag- 
CO ■ netic arcade stretching and plasmoid formation induced by the differential rotation of the solar 
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CO photosphere or of an accretion disk. 
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We present a new class of soliton-like solutions of the Grad-Shafranov equation. The 
Grad-Shafranov (GS) equation governs the equilibrium conditions of a magnetic field em- 
bedded in a plasma. 

Countless applications rely on the GS equation to determine the equilibrium steady state 
condition jlj|. The equilibrium field in magnetic fusion devices is computed with the GS 
equation. In the present letter, the attention focuses more on problems where the GS 
equation can be of importance for space and astrophysical systems. 

In the solar corona, field lines emerging from the photosphere form arcades whose evolu- 

Q 

tion is crucial in solar flares and coronal mass ejections [2j. Similarly, magnetic fields tied 
to an accretion disk rotating around a massive object (such as a star or or a black hole) can 
form magnetic arcades j^J . In both cases, photospheric motions or the differential rotation of 
the accretion disk cause a shear of the footpoints of the magnetic arcade. The shear motion, 
in turn, injects toroidal magnetic flux and causes the expansion of the arcade away from the 
photosphere or from the accretion disk. 

In this scenario, the role of the GS equation has been traditionally tied to the under- 
standing of the expansion of the magnetic arcade. A series of papers have investigated the 
expansion phase using self-similar solutions of the GS equation, proving that the field ex- 
pands away from a differentially rotating surface at a characteristic angle of 60 degrees 
However, the applicability of the GS equation to this scenario is intrinsically limited by 
two factors. First, the GS neglects inertia effects that can be important in determining the 
expansion of the arcade. Secondly, the presence of dissipation processes (such as anomalous 
resistivity or kinetic effects) can lead to reconnection processes that alter the topology of the 
expanding arcade. Several papers have investigated this process proving that reconnection 
can lead to the formation and ejection of a plasmoid, i.e. a blob of plasma encircled by field 
lines at least partially detached from the original footpoints. In the solar case such plasmoid 
can be interpreted coronal mass ejection. 

Recently it has been pointed out that the process of arcade expansion and plasmoid 

n 

ejection can repeat itself in time at the same location [5|. There is observational evidence 
to support this claim. In the solar case, it has been observed that flares can repeat at the 
same location on a cycle of several hours to few days In the case of accretion disks, 
the repetitive ejections of plasma "bullets" has been observed for SS433 p and jets emitted 
from accretion disks often present a knotty structure that suggests the presence of multiple 
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islands [7]. 

In the present letter, we discuss a new model of such ripetitive series of plasmoids. The 
model ineterprets such observational occurencies in the mathematical framework of soliton 
theory. A jet of plasma composed by a series of plasmoids is represented mathematically 
as a soliton being propagated from the source (the photosphere in the solar case and the 
accretion disk for the astrophysical case). 

Below, we show that a mathematical analogy exists between the GS equation and the 
cubic Schrodinger (CS) equation typical of the soliton theory Such partial analogy is 
used to find a class of exact solutions of the GS equation that have the same mathematical 
structure of the soliton solutions of the CS equation. Finally we discuss the likelihood of 
such solutions existing in practice. 

Grad-Shafranov Equation: Magnetic equilibria are computed using the GS equation. The 
GS equation is derived straight-forwardly from the momentum equation assuming equilib- 
rium conditions and neglecting inertia: 

J x B = Vp (1) 

Assuming a Cartesian geometry with axis (x,y,z), we assume all quantities to depend 
only on (x, y) and assume the third dimension to be a symmetry direction. In most of the 
applications, axisymmetric coordinates would be a more appropriate choice, but not within 
the scope of the present work. The use of the axisymmeric coordinates would obscure 
the mathematical equivalence that we want to prove between the GS equation and the CS 
equation. The extension of the derivations below to axisymmetric coordinates will be the 
topic of a future paper. 

With the geometry chosen above, a solution of eq. (J^) can be found in the form jjj: 

B = z x W + B z z (2) 

where z is the unit vector in the ignorable direction z. From the magnetic field given in 
eq. (j2J), the current follows immediately: 

J = V x B = VB Z x z + V 2 #z (3) 

Substituting the expression for the current and the field, the balance eq. (JTJ) assumes a 
new form: 

{VB Z x z) • Wz = + V 2 tf W + B Z VB Z (4) 
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From this rather involved equation, the GS equation follows, noting that the left hand 
side is directed along the z direction and the right hand side is directed on the (x, y) plane. 
It follows that the two sides of the equations must both be zero for the complete equation 
to be satisfied: 

(W x VB Z ) • z = 

(5) 

Vp + V 2 #W + B Z VB Z = 
where the left hand side of eq. (j3J) has been manipulated into a more convenient form. 

The two equations must be both satisfied to obtain a complete equilibrium. Note that 
the first equation simply states that the gradients of B z and of $ must be parallel. This 
requirement is typically enforced by requiring that B z is itself a function of \1/ (the same 
being also true for p) : 

'p = </(*), B z = f(V) 

(6) 

V 2 M> = -g' - ff 

v 

where the prime represent derivative with respect to ^. Equation © is the classic form of 
the GS equation. 

Typically it is convenient to work on the classic form © of the GS equation obtained 
above, but we will show that it has a limitation: it does not allow to consider the extension 
of \1/ in the complex plane. Often in mathematical physics, it can be advantageous to extend 
an equation in the complex plane to obtain its solutions, returning to the real axis after 
the solution is obtained. To obtain the class of soliton-like solutions derived below, that is 
indeed the procedure to follow. To that end, the original form (J5J) of the GS equations is 
to be preferred, as the independent variables (x, y) remain real even when the dependent 
variable \1/ is extended in the complex plane. In the classic form of the GS equation, the 
extension of \l/ to the complex plane would entice the consideration of complex functions of 
complex variable. The use of the original form (JHJ) requires only the consideration of complex 
functions of real variables. 

Analogies for the Grad-Shafranov Equation: When studying the GS equation, one is 
confronted with two possible analogies. One is well known and has been often used in 
previous works: the analogy with the Helmholtz problem of mathematical physics. The 
other is new and will be presented here for the first time: the analogy with the CS equation. 
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Eq. (jUJ) when considered in the traditional way as an equation for a real function of a 
real variable ^(x, y) is an application of the Helmholtz problem: 

AT = AT (7) 

for the eigenvalue A and the eigenf unction T. 

In studies of the GS equation applied to space and astrophysics systems, the analogy with 
the Helmholtz problem (JZj) has been very fruitful in studying arcade expansion jj]. However, 
all solutions provided by the analogy with the Helmholtz problem do not include the presence 
of reconnection. As noted above, reconnection is believed to be present as suggetsed by 
simulation and supported by observation. To arrive to a model that incorporates the effects 
of reconneciton and the presence of plasmoids, we need to abandon the analogy with the 
Helmhotlz equation and propose a new one: the analogy with the CS equation. 

The inductive process that leads from the Helmholtz problem to the CS equation is 
typical of the discipline of nonlinear propagation in optical waveveguides. Wave propagation 
in optical waveguides is governed by a Helmholtz equation for the electromagnetic field [ll]] : 

= -k 2 n 2 ^ (8) 

where n is the refractive index of the media, k the wavenumber and \1/ is a complex function, 
interpreted physically as a polarization component of the electromagnetic field. In non- 
homogeneous media n is a function of space and in non-linear optical materials it is a 
function of the electric field amplitude: n — n(\^f\). In the most common case of the Kerr 



effect jl(3], the dependence of the refractive index upon the amplitude of the field is quadratic: 
n = riQ + n2|\l/| 2 . Following the classic textbook derivation ijj, we assume a 2D system 
(in analogy with the system under consideration for the GS equation) and we rewrite the 
unknown function using a leading harmonic dependence of the field along the propagating 
direction (that we assume to be y in the present case): 

*(x,y)=<f(x,y)e- jknoy (9) 

where j is the imaginary unit. Absolute generality is still provided by retaining the most 
general dependence of the factor ip(x,y). For the new unknown function ip, the complex 
Helmholtz equation (JBJ) yields: 

|V + 0-2^1^-^ = (10) 



Traditionally, a further approximation is made ^(|, assuming (n 2 — 71q) = 2n (n — n ) 
and requiring the applicability of the paraxial (or Fresnel) approximation that prescribes 
| <C 2kn \^\. The final equation typically used to study light propagation in non linear 
media becomes: 

none others than the cubic Schrodinger (CS) equation A well known theoretical result 



is that the CS equation above admits a soliton solution of the form 
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* = ^ sech © exp (^§) (12) 

where y Q is a free parameter and |y?p| 2 = 2/n n 2 y 2 . 

The formal analogy of the original Helmholtz equation in the complex plane used to study 
light propagation with the GS equation allows to derive a new class of soliton-like solution 
similar to eq. (refsolitonCS). 

Soliton solutions of the Grad-Shafranov Equation: To use the analogy with the CS equa- 
tion outlined above, we consider the GS equation in its original form (pj) and assume the 
following choice of the free functions B z and p: 



B Z VB Z = afoVV 



Vp = a£|#| 2 *V* 



(13) 



The GS equation (J5J) is extended in the complex plane and complex functions ^ of real 
variables (x, y) are sought as solutions. The choice above (fTB^) is a possible and legitimate 
choice of the free functions in the GS equation. Although at first sight it might appear 
arbitrary, the choice selected above has a simple physical meaning: it ensures that VB 2 /Vp 
is independent of the y coordinate along which the soliton-like solution propagates. With 
this choice the first of the GS eq. (0) is automatically satisfied and the second becomes 
identical to eq.fjHJ) with a non-linear refractive index of the form prescribed by the Kerr 
effect: 

d 2 V d 2 ^ 2/1 IT|2 . 

= + (14) 

A soliton solution can be readily obtained in the form: 

V(x,y) = y p sech{x/L)e- j{ao+1/2aoy ° )y (15) 
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FIG. 1: Contour plot of the real part of the flux function, corresponding to the section of the 
magnetic surfaces on the (x, y) plane. 

with L = {l/yl + l/4a^)- 1/2 and # p = V2/a L. 

Note that the soliton solution (|15jl is exact and no paraxial approximation has been made 
in the present case, since the exact solution can be found even for the exact GS equation 
and the physical significance of the paraxial approximation is peculiar to optical waveguide 
propagation, and of no meaning in the case of the GS equation. 

Figure H shows the contour lines of the real part of the flux function physically such 
contour lines correspond to the section of the magnetic surfaces on the (x, y)plane. The field 
lines belong on the flux surfaces but have a third component away from the plane and given 
by B z : 

B z = a Q m (16) 

The magnetic field topology presents a series of magnetic islands piled one above the 
other. Physically, this solution can be considered the end result of a repetitive series of 
processes of arcade stretching and plasmoid ejections that compose a jet of plasma and 
magnetic field emitted from a differentially rotating surface (such as the photosphere or an 
accretion disk as described above). 

It is interesting to observe that the primary balance of force in the present case is between 
the magnetic tension of the field lines and the plasma pressure. Using the solution (|15p. the 
initial choice (fT3]l for the free functions in the GS equation yield an equation for the pressure 
perturbation: 

A5p = V ■ (a 2 \ty\ 2 ^W) (17) 



FIG. 2: Dependence on x of the pressure (solid) and of the flux function (dashed), ao = 1, yo = 1- 
In astrophysical systems, it is reasonable to assume a background pressure that is being 

I! rrm 

perturbed by the magnetic structure |3|]. Eq. (jTTj) can be readily solved observing that the 
y dependence of 5p is harmonic: 



Substituting into eq. (|T7|) . a simple ordinary equation is obtained for the x-dependent factor 
X of the pressure perturbation: 



that is solved numerically. 

Figure 121 shows the x dependence of the pressure \. The y dependence is harmonic and 
is given by eq. (fTSj) 

Existence in Nature of soliton-like solutions: The solution derived above is an exact and 
legitimate solution of the GS equation. It has been derived in the complex plane and to 
obtain a physically viable solution one has simply to take either the real or the imaginary 
part. 

Once we have shown the mathematical existence of soliton-like solutions, we need to 
consider the possibility that in practice they can be identified as realistic descriptions of 
processes observed in nature or generated in laboratory. Here we propose that soliton- 
like solutions can explain the observed structure of astrophysics jets. Astrophysics jets are 
observed to be emitted from rotating accretion disks. The dynamo effect inside the accretion 
disks generate a magnetic field that is then emitted with the flow of the jet |2|. Observations 
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show that often jets present a so called "knotty" structure where the jet is composed by a 
number of islands (knots) 7]. This feature suggests the possibility to model the magnetic 
structure of jets using multiple islands solutions of the Grad-Shafranov equation. While 
the solution-like solutions are by no means the only solutions presenting multiple islands 
(e.g. for another example), they present the right qualitative behaviour observed in 

astrophysics jets. Two issues need to be considered to propose credibly the applicability of 
soliton-like solutions. 

The first question is whether a soliton-like solution is stable. Previous studies have 
considered the stability of equilibria composed by periodic magnetic islands, proposing the 
existence of the coalescence instability 12|. For sufficiently large aspect ratio islands, i.e. for 
islands stretched along the axis of the jet (the y axis of the solution above) the instability 
is greatly reduced allowing the soliton-like solution to remain stable for long times (that 
observationally turn into long distances, recalling that the jets are moving at relativistic 
speeds). One can indeed propose that the eventual onset of this instability can explain why 
the jets after remaining collimated for kiloparsec distances eventually destabilize in large 
magnetic clouds. 

The second question is whether a soliton-like solution can be created by a disk. To answer 

this question one should do a simulation of the accretion disk and of the jets emitted from 

them. A complete treatment is extremely challenging and beyond the scope of the present 

work (if at all possible). Yet limited studies that focus on just a part of the problem are 

possible and have been performed. The creation of repetitive islands similar to the islands 

featured by the soliton-like solution have indeed been observed in simulations of a problem 

1 

closely related to jet formation: the study of plasma emitted from the solar corona p. And 
indeed the ejection of multiple magnetic islands has been observed for the jet from SS433 [f| 
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